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Abstract 



We analyze the non-Gaussian primordial fluctuations which are inescapably contributed 
by scalar fields $ with vanishing expectation values, (<fr) = 0, present during inflation 
in addition to the inflaton Geld. For simplicity we take $ to be non-interacting and 
minimally coupled to gravity. $ is a Gaussian variable, but the energy density fluc- 
tuations contributed by such a field are x 2 -distributed. We compute the three-point 
function $,J; S ^ (£) for the configuration of an equilateral triangle (with side length i ) and 
the skewness (S^), i.e. the third moment of the one-point probability distribution of 
the spatially smeared energy density contrast 5r, where R is the smearing scale. The 
relative magnitudes of the non-Gaussian effects, [^ iV ^] 1/ ' iV /[^ 2 -'] 1 ^ 2 , do not grow in time. 
They are given by numerical constants of order unity, independent of the scale I. The 
" bi-skewness" (5^5 s ) is positive. For smearing lengths R <C S this shows that in our 
model (in contrast to Gaussian models) voids are more quiet than high-density regions. 

PACS numbers: 04.62. +v, 98.80.Cq, 98.65. Dx 
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1. Introduction and Conclusions 

Gaussian primordial fluctuations are predicted by standard inflationary models [1] . 
The data on the anisotropies of the cosmic microwave background taken by the COBE 
satellite [2] are not sufficient to give a significant test for this prediction of a Gaussian 
distribution [3]. Large non-Gaussian effects are observed in the large-scale structure 
of the galaxy distribution, e.g. in the one-point probability distribution of the matter 
density p, smeared over a window scale R, where R is in the range 5 to 20 h~ x Mpc 
[4]. One inevitable source for this non-Gaussianity is the non-linear character of the 
gravitational evolution, irrespective if the initial fluctuations are Gaussian or not. 

To test the hypothesis of Gaussian primordial fluctuations, one needs one or several 
specific non-Gaussian models. One class of models which give rise to non-Gaussian 
primordial fluctuations involves cosmo logical defects, like cosmic strings or textures [5]. 

In this paper we discuss a different origin of non-Gaussianity in the primordial 
fluctations, and we analyze a minimal non-Gaussian model, which has been introduced 
in ref. [6]. Non-Gaussian primordial fluctuations are inescapably contributed by scalar 
fields $ with vanishing expectation value, (<fr) = 0, which are present during inflation in 
addition to the inflaton field. For simplicity we consider $ to be a non-interacting scalar 
field, minimally coupled to gravity. $ is a Gaussian variable. Since the energy density p 
of this non-inflaton field is bilinear in the fluctuating field p is x 2 -distributed. (Recall 
that the x 2 -distribution is defined as the probability distribution of the mean square 
error of a Gaussian variable). In contrast the energy density fluctuations of the inflaton 
arise from the interference terms between the background field and the fluctuating part 
of the inflaton field. Therefore the energy density fluctuations are linear in the inflaton 
field fluctuations, and the inflaton energy density fluctuations are Gaussian. 

The order of magnitude of the fluctuations at the second horizon crossing con- 
tributed by non-inflaton massless scalar fields is Sp/p ~ Hf /Mpj, where Hi is the Hub- 
ble constant during inflation. For massive, unstable scalar fields with masses m$ <C Hi, 
5p/p acquires an enhancement factor ^/m$/r$, where r$ is the decay rate [7]. In 
this short note we consider only the massless case, in refs. [6,7] we have analyzed both 
m$ = and m$ 7^ 0. In previous analyses of non-Gaussian fluctuations in inflationary 
models with two scalar fields [8] the physical origin and the resulting character of the 
non-Gaussianity are very different from the ones presented here. 

Observable measures of non-Gaussianity are the expectation values of products 
of N energy density contrasts S(x) for N > 2. In this paper we restrict ourselves 
to N = 3. We first (in sect. 2) compute the three-point function £( 3 )(a?i, x<i, £3) = 
(S(xi)d(x2)S(xs)) for the configuration of an equilateral triangle with sides of length I, 
i.e. CA S H^)- As a second observable we compute (in sect. 3) the skewness (8 R ), i.e. the 
third moment of the one-point probability distribution of the spatially smeared energy 
density contrast 5r := J d 3 x 8{x)Wr{x) with a window function Wr{x) of smearing 
scale R. As a third measure of non-Gaussianity we consider (in sect. 4) the "bi- 
skewness" (8 R 8 S ) in the case R <C S. 
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The non-Gaussianity of the primordial fluctuations in our model has a simple 
structure. For a x 2 -distributed variable, such as the energy density contrast 5, the 
expectation value of a product of N 5's is of the order of the (N/2) th power of the 
expectation value of two <Ts, e.g. for the iV-point function we find a of the order of 

[£(2)jiV/2_ Thig 

is so because any £( N ) must be built out of products of N/2 two-point 

functions of 

The relative magnitude of the non-Gaussian effects in £( N ^ is given by 
|-£(7V)ji/A^|-£(2)ji/2_ rj ur j n g the cosmological time evolution the perturbations <^ 2 ) grow 
for fixed comoving £. But in our class of models (with its linear dynamics) the ratios 
|-£(JV)ji/.Ny|-£(2)ji/2 are independent of time. They are given by numerical constants of 
order unity. 

The cosmological model considered is an inflationary de Sitter era, followed by a 
radiation dominated era. The quantum state of the field $ is the Bunch-Davies state 
initially, during inflation. The scalar field $ (not the inflaton) evolves in this background 
curved space-time, i.e. the back reaction of $ on the geometry is neglected. Therefore 
gauge ambiguities are eliminated [9] . Since we do not treat the decay of the scalar field 
(into radiation or other forms of matter), we make our predictions for the end of the 
radiation dominated era. 

At the end of the radiation era (the time of matter and radiation equality) the 
model contains two characteristic scales, the Hubble parameters H\ and -Heq- For 
cosmologically relevant length scales the quantities £a^C0 an d each contain two 
domains according to whether I or R are smaller or bigger than -Heq- On sub-horizon 
scales, H^q£ <C 1, the three-point function increases logarithmically towards the smaller 
scales, £ A (3) W ~ |logH£| 3 . On super— horizon scales, H-^qI 3> 1, we obtain £a (-0 ~ 
(Hi)- 6 . We find the scaling laws £ (3) /[£ (2) ] 3/2 = 5/VG for sub-horizon and = ^8/27 for 
super-horizon scales. We obtain (8 R ) / (8 R ) 3 ^ 2 = 10/(67re 3 ) 1 / 4 for the ratio of skewness 
and variance of 5r with smearing scales HeqR <C 1. The " bi-skewness" (S R 8 S ) is 
positive for all scales R and S. For R <^ S this shows that in our model the small-scale 
variance (8 R ) observed in a sample within a void of size S ^> R is expected to be smaller 
than the small-scale variance observed in a high-density region of size S. 

The non-linear gravitational evolution of Gaussian primordial fluctuations produces 
another simple relation between expectation values of iV 5's in lowest non-vanishing 
order perturbation theory [10]. In this approach the iV-point functions are of 
the order [^ 2 - ) ] Ar_1 , completely different from the relation valid for a x 2 -distributed 
quantity. 

According to refs. [11,12] the observed iV-point correlation functions of galaxy 
counts for iV = 3, 4 in the range 1 to 50 h' 1 Mpc are consistent with non-linear 
gravitational evolution from Gaussian primordial fluctuations. In order to see how 
much room the observations leave for non-Gaussian primordial fluctuations, it would 
be necessary to apply the non-linear gravitational evolution to these fluctuations. 
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2. The Three— Point Correlation Function 



The background geometry is a Friedmann-Robertson- Walker space-time with spa- 
tially flat sections, p to t = Pcrit- With conformal time r\ we have ds 2 = a(rj) 2 {drf — 
dxidx 1 ). The Hubble constant during inflation is denoted by H\. The transition be- 
tween inflation and the radiation era can be approximated as instantaneous, because the 
physical transition time is much shorter than the characteristic time for the evolution of 
the cosmologically relevant modes, which have A p h ys > e 70 ^ -1 . We fix the coordinates 
in such a way that at the time r]i of the transition {i.e. at the end of inflation) we have 
r]i = Hi~ x and a(r]i) = 1. This fixes the scale factor as 

air]) = j- — . (r]<H7 1 ) 
yu (2-H lV ) w " 1 ' (2.1) 

a( V ) = Htn (V > Hr 1 )- 

The Hubble parameter H(rf) during the radiation era reads H = a~ 2 H\. 

We consider a non-interacting, massless, neutral scalar field $ with minimal cou- 
pling to gravity. Its energy density p is 



P = 



2a 2 

The energy density contrast contributed by the scalar field considered is 



(2.2) 



S(x) = J-(p( x )-(p», (2.3) 

Ptot 

where p tot = 3H 2 M^/8n. 

The equal time three-point function of the energy density contrast is 

£ (3) (a-i,f 2 ,£ 3 ) = (J2|(y(xi)<J(f 2 )<y(f 3 )|fi). (2.4) 

The quantum state |0) is initially (during inflation) the Bunch-Davies state: Every 
observationally relevant mode, which has today (A p h ys )o < Hq 1 , had at early times 
in inflation -R//c p h ys 2 - ► and is taken to be initially in the Minkowski vacuum state. 
Since the quantum state \Q) is translation invariant, depends only on the mutual 
physical separations of the three points Xi. 

We work with a normal ordered energy density operator N[p], and we use the fact 
that any normal ordering N gives the same result for the correlation functions. This is 
so because the difference of two normal orderings is a c-number, which drops out in the 
density contrast, eq. (2.3), and a fortiori in the correlation functions. 

Since the results are independent of the normal ordering, we can make the most 
convenient choice, the one adapted to the quantum state For the mode expansion 
of the field operator $, 

/d s k 
J2^)3 l a kn^^) + (2.5) 
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we choose the expansion in terms of those evolving basis modes p^iv, %) which at very 
early times, when A p h ys <C H^ 1 , behaved as ~ ex.p(ikx — ikrf). It follows that 



ar\Q) = 0, Vfc. 



(2.6) 



One could write more explicitely (p^iVi x) = (p^ in \rj, x), = a^ in \ and \Q) = \BD, in) 
(where BD stands for Bunch-Davies). For more details about the physics and the 
functional form of the evolved modes we refer to refs. [6,7]. 

In [6] it was shown that the Wightman function W(x, x') of the state 



/d 3 k 
^3 n {x) n {x'T, 



(2.7) 



is extremely useful, because the equal time two-point function of the energy density 
contrast ^ 2 \x, x') = (Q\8(x)5(x')\Q) in the state \Cl) can be expressed in terms of 
derivatives of W(x, x') as 



(2a 2 p t ot) ai ,a 2 =0 



[W(x 1 ,x 2 ) 

,a,\a,2 J 



(2.8) 



The derivatives with index i act exclusively on the coordinates of the point Xi. For the 
three-point function ^ 3 ~> the analogous expression reads 

^ S \xi,X2,X 3 ) = 

8 3 

-3 Y W(x U X 2 ), aia2 ■ W(x 2 ,Xs), a2 a 3 ■ W(x3,Xi) itt3ai 

(2a 2 p to t) ai) a 3 ,a 3 =0 Vl=V2=V3 

(2.9) 

Cosmologically relevant separations of the three points, today 1 to 3000 h~ x Mpc, 
correspond to physical point separations at the end of inflation much larger than the 
Hubble radius H^ 1 at that time. Therefore we can use an approximate Wightman 
function for the state \CL) during the radiation era [7], 



W(x,x') = H * 



1 



(2tt) 2 24r]7]'A 



A 3 _ + log|^iA_+| + A+ _ log \HiA + . 



- A 3 __ logl^A..) - A3 + log|H I A ++ | 



(2.10) 



The comoving distance of the two points is denoted by A~ := \x — x'\, and we have 
introduced the notation A±± : = A" ± rj ± rj'. 

The three points x% in eq. (2.4) define the corners of a triangle. We consider 
the special case where the three points form an equilateral triangle, and we denote 
the corresponding three-point correlation function by £a^C0> where £ is the physical 
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distance between any two of the points. The model contains two intrinsic length scales, 
namely H^ 1 and H~ x . For cosmologically relevant point separations, i 3> aH^ 1 , the 
three-point function contains two domains, according to whether i is bigger or smaller 
than H~ x . For sub-horizon length scales, Hi <Cl, we find 



On super-horizon scales, Hi 3> 1, we obtain 



Hi 
Mpi 



, Hi 
log — 



(2.11) 



Ca (3) W * ( 



( 2 



Hi 



V37ry V m pi/ CW 



(2.12) 



Comparing eq. (2.8) for £( 2 ) with eq. (2.9) for £( 3 ) we expect that there should be 
a scaling law £ A (3) (0 = <D[^ 2 \i)] 3 / 2 . Using the results of ref. [7] for £ (2) (£) we find 



^(2)(£)]3/2 



-= for £ < 7i _1 
v/6 



for 




(2.13) 



On sub-horizon (7Y£ <C 1) and super-horizon scales (7i£ ^> 1) we have scaling laws. 
The constant coefficients are different in the two regimes, such that in the vicinity of 
i = H~ x the ratio £a (3) (0/[£ (2) W] 3/2 is slightly ^-dependent, see figure. 



3. The Skewness 

We introduce the spatially smeared energy density contrast 

S R := J d s x S(x)W R (x). (3.1) 
For convenience, we take a Gaussian window function W R of scale R, 

We compute the second and third moments, (5 R ) and (8 R ), of the one-point probability 
distribution for 5 R . Note that one-point probability distributions inescapably need a 
smearing scale R, both in the analysis of observed galaxies and in quantum field theoretic 
computations. The reason is that the observed mass density in galaxy counts is the sum 
of Dirac delta functions (galaxies are treated as points), which cannot be taken to the 
N th power without smearing beforehand. In quantum field theory, on the other hand, 
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the universal short-distance behaviour, e.g. i^ 2 \£) ~ £~ 8 , necessitates smearing before 
squaring. The cosmologically relevant term has a short-distance behaviour ~ (log7i£) 2 , 
which also necessitates smearing. 

The variance (8 R ) = J d 3 x J d 3 yW R (x)W R {y)^ 2 \x,y) is given by a double in- 
tegral. Since the equal time two-point function £( 2 ) can depend only on the physical 
separation £ of the two points, one integration is trivial. We obtain 

The skewness (5 3 R ) = J d 3 x J d 3 y J d 3 zW R (x)W R (y)W R (z) ^ 3 \x,y, z) is given as a 
triple integral. The equal time three-point function ^ 3 \ which depends only on the 
mutual physical separations of x, y, and z, can be characterized by any three quantities 
defining the triangle given by the points. We choose two sides of length £ and £' and 
the angle 9 in between, ^ 3 ^ (£,£'; cos 9). One integration is again trivial, and we arrive 
at 

1 f°° dl f°° d£' f 1 n q £ 2 +£' 2 -£f cose 

®=V^i Tj ^jj^^e — {«(^;cc..). 

(3.4) 

Let us first concentrate on the case of sub-horizon smearing scales, HR <C 1. For 
sub-horizon separations, £,£' <C 7Y _1 , the three-point correlation function varies loga- 
rithmically, ^ 3 \£J';cos9) ~ ^if I 6 Mp 1 6 (2/37r) 3 |log^logWlogW|, where 

£" = \£ — £'\. The integral (3.4) therefore is strongly peaked at values £ = £' = 
— cos 9, and we obtain 

- 5 S {I) 3 {wS llosHR ^ (HR<<1) - (3 ' 5) 

Comparing the variance [7] and skewness we obtain the numerical constant in the general 
scaling law discussed in the introduction, 



10 (HR«1). (3.6) 



[(5 2 R )} 3 / 2 ~ V 67re 3 

For super-horizon smearing scales, HR ^> 1, the exponential in eq. (3.4) can be 
dropped, and the skewness (5 R ) is approximately 

^- 8 i£ c >{0^ <«*» ^ 

where the numerical constant C3 is defined by the integral of ^ 3 \ 

/ a / a ' 5<3,( "'°-(t) 3 (^) 6 §' < 37b > 
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For HR ^> 1 the general scaling law is violated. The reason for this is simple: If the 
physical separations of the points involved in ^ 2 ' 3 ^ are larger than H~ l , both correlation 
functions show a steep fall-off, ^ 2 \£) ~ (ftf)" 4 and ^ 3 \£,£' ; cos6>) ~ {H 3 ££' £")- 2 . The 
integrals (3.3) and (3.4) are therefore strongly peaked at values for £, £', and I" of the 
order of 7i _1 , independent of the scale R ^> l~L~ l . But the /^-dependence from the 
prefactors in eqs. (3.3) and (3.4) gives a scaling behaviour such that another ratio, 
(5 R )/[(5 R )] 2 , is independent of R. This reflects the fact that (5 R ) involves three non- 
trivial integrations, while for (8 R ) six are needed. 

The fact that the skewness is positive in our model (with its linear dynamics) 
is totally independent of the kinematical effect, p > 0, 5 > — 1, which forces large 
fluctuations (needing non-linear dynamics) to be associated with positive skewness. 



4. The Bi— Skewness 



In the analysis of expectation values of three observables in position space, we now 
consider, after %A 3 \£) and (S R ), the " bi-skewness" (S R S S ), where the two smearing 
scales R and S (for density contrasts at the same point) are different. This measure 
of non-Gaussian effects appears to be both robust (observationally) and discriminating 
(with respect to models) . For scales R <C S the bi-skewness (S R S s ) provides the answer 
to the following physical question: Is the small-scale variance (5 R ) observed in a sample 
within a void of size S (S 3> R) expected to be smaller, equal, or larger than the small- 
scale variance observed within an extended high density region of size S ? For the cosmic 
microwave background anisotropies the physical question is: Are the small-angle (R) 
anisotropies observed within a moderately large patch (S) over the South Pole expected 
to be the same as those observed within a moderately large patch (S) over the North 
Pole, in view of the fact that the average temperatures over the two different samples 
are different ? For Gaussian fluctuations the small-scale variance is expected to be the 
same for samples within large voids and for samples within large high-density regions. 
For our model of non-Gaussianity (with its linear dynamics), the voids are more quiet 
than the high-density regions (with respect to the small-scale variance). For R <C S 
we find ($%S S ) > 0. We obtain the integral expression in the limit R <^ S 

1 r°° d£ f°° d£' f 1 a q P i' 2 -W cose 

^ R S/ VSnR^S 3 Jo £ Jo £' J-i (4.1) 

x^ 3 )(£,f;cos6>). 

We insert the Wightman function (2.10), and for sub-horizon smearing scales, 
R < S < 7Y" 1 , we obtain 

The physics reason why the small-scale (R) variance is smaller in samples (of size 
S ^> R) from voids than in samples from overdense regions is easily seen, when we 
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consider only two modes for simplicity, one mode with ~ R, the other mode with 
of the order of the separation of the two samples. Since dp contains interference 
effects between the two modes, the void is less noisy than the overdense region. 

In the opposite limit, when the two scales obey R ^> S, the bi-skewness (S^dg) 
is an uninteresting quantity, since it contains no new information in addition to the 
skewness (<5#). 
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